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Abstract. Corrections to the relativistic theory of orbits are discussed considering higher 
order approximations induced by gravitomagnetic effects. Beside the standard periastron 
effect of General Relativity (GR), a new nutation effect was found due to the c"^ orbital 
correction. According to the presence of that new nutation effect we studied the gravita- 
tional waveforms emitted through the capture in a gravitational field of a massive black 
hole (MBH) of a compact object (neutron star (NS) or BH) via the quadrupole approxima- 
tion. We made a numerical study to obtain the emitted gravitational wave (GW) amplitudes. 
We conclude that the effects we studied could be of interest for the future space laser inter- 
ferometric GW antenna LISA. 

Key words, theory of orbits - gravitomagnetic effects - stability theory - gravitational 
waves. 



1. Introduction 

The magnetic field is produced by the motion 
of electric-charge, i.e. the electric cuiTent. The 
analogy with gravity consists in the fact that a 
mass-current can produce a "gravitomagnetic" 
field. The formal analogy between electromag- 
netic and gravitatio nal fields was explored by 
IA. Einstein! (1 191 3h. in the fram ework of GR, 
and then by'H.Thirrins' ('1918V It was shown, 
by [j. Lense and H.Thirring ( 1918), that a ro- 
tating mass generates a gravitomagnetic field, 
which in turn, causes the precession of plan- 
etary orbits. We want to study how the rela- 
tivistic theory of orbits and the production of 
GW is affected by gravitomagnetic corrections. 
The coiTections, which are ofF-diagonal terms 
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in the metric, can be seen as further powers in 
the expansion in c"' (up to c^^). Nevertheless, 
the effects on the orbit behavior involve not 
only the precession at peri-astron but also nu- 
tation corrections. Our approach suggests that, 
in the weak field approximation, when con- 
sidering higher order corrections in the equa- 
tions of motion, the gravitomagnetic effects 
can be particularly significant. In systems ap- 
proaching to strong field regimes, these cor- 
rections give rise to chaotic behaviors in the 
transients dividing stable fr om unstable orbits 
Is. C apozzie llo et aTl (12009 ^. In general, such 
contributions are discarded since they are as- 
sumed to be too small but they have to be taken 
into account as soon as the ^ ratio is significant. 
It is possible to take into account two types of 
mass-cuiTent in gravity. The former is induced 
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by matter source rotations around the center of 
mass: it generates the intrinsic gravitomagnetic 
field which is closely related to the angular mo- 
mentum (spin) of a rotating body. The latter is 
due to the translational motion of sources. 



2. Gravitomagnetic corrections 

Starting from the Einstein field equations in 
the weak field approximation, one obtains the 
gravitomagnetic e quations and then the correc - 
tions in the metriclcW. Misner et al. I (Il973l) : 
Is. Capozziello et al. I ( l200'll) : 



, / 20\ , , 85/, y' 
ds-^ = ll + — \ c^dr —cdtdx-' 



( 2<D\ 

1 — I dijdx'dx^ . 



(1) 



By calculating the affine connection related to 
the metric ([T]i, one also obtain the geodesic 
equations 



x" + r^.x^x" = , 



where the dot indicate differentiation with re- 
spect to the afiine parameter In order to put in 
evidence the gravitomagnetic contributions, let 
us explicitly calculate the Christoffel symbols 
at lower orders. By some straightforward cal- 
culations, one gets 
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In the approximation which we are go- 
ing to consider, we are retaining terms up 
to the orders O/c^ and jc'. It is im- 
portant to point out that we are discard- 
ing terms Hke ((D/c4)5<l)/5x*, (y Vc-'*)5<D/5x*, 
('^lc^)dV''ldxK iV''/c^)dVj/dx' and of higher 
orders. Our aim is to show that, in several cases 
like in tight binary stars, it is not correct to 
discard higher order terms in v/c since phys- 
ically interesting effects could come out. The 



geodesic equations up to c corrections are 



then 
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for the time component, and 
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for the spatial components. Considering only 
the spatial components, we obtain the orbit 
equations. Calling t/Ze„c/,v/ - dijdxjdxj, ande*^ - 

, I- 

we have, in vector form. 



dl, 



de 



(2) dleucUd 



= -4[V<l)-e(e-V<l))]-i- 



^ [e A (V A V)] 
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The gravitomagnetic term is the second one in 
Eq.© and it is usually discarded since consid- 
ered not relevant. This is not true if v/c is quite 
large as in the cases of tight binary systems 
or point masses approaching to BH. From the 
above equations we can write the Lagrangian 
and derive the orbital equations of motions 
start ing from the Euler-Lagran ge equations 
(see Is. Capozziello et alTI (12009 ). Our aim is 
to study how gravitomagnetic effects modify 
the orbital shapes and what are the parame- 
ters determining the stability of the problem. 
The energy, the mass and the angular mo- 
mentum, essentially, determine the stability. 
Beside the standard periastron precession of 
GR, a nutation effect is induced by gravito- 
magnetism and stability depends on it. The so- 
lution of the system of differential equations 
(ODE) of motion presents some difficulties 
since the equations are stiff. For our purposes, 
we have found solutions by using the so called 
Stiffness Switching Method to provide an au- 
tomatic mean of switching between a non-stiff 
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and a stiff solver coupled with a more con- 
ventional explicit Runge-Kutta method for the 
non-stiff part of differential equations. Time se- 
ries of both r{t) and r{t) together with the phase 
portrait r(t), r{t) are shown assuming as initial 
values of the angular precession and nutation 
velocities integer ratios with the radial veloc- 
ity. In Fig. [U we show as example the one with: 
^ = and 61 = ^(p 



M=1 En=0.95 r0=20 -=— -= — 
i 10 r 10 




Fig. 1. Plots along the panel lines of: r{t) (upper 
left), phase portrait of r{t) versus r{t) (bottom left), 
r(t) (upper right) and r{t) (bottom right) for a star of 
IMq. The examples we are showing were obtained 
solving the system for the following parameters and 
initial conditions: fi » IM©, E = 0.95,(^o = 0, 
Oq = f , 6o = -j^^o, 00 = -j^ro and tq = --^^ and 
To = 20fi. The stiffness is evident from the trend of 
r(t) and r(t) 



En= 0.95 r„= 20 r = -0.1 tb = 0.01 6 = 0.001 
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Fig. 2. Plots of the basic orbits with the associated 
field velocities in false colours 

Orbil with Gravitomagnetic effect(Red] and without Gravitomagnetic(black} 
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2.1. Orbits with gravitomagnetic 
corrections 

In this section we show some examples of 
orbits. In Fig|2] are plotted some basic orbits 
with the associated field velocities in false 
colours. Then in Fig.[3]we show the orbits with 
gravitomagnetic correction (red-line) and with- 
out gravitomagnetic correction (black-line). 
Finally in the Fig. |4] there is the phase portrait 
with (red-line) and without (blue-line) gravit- 
omagnetic orbital correction respectively. The 
example we are showing was obtained solving 
the system for the following parameters and 
initial conditions: fi « IM©, E = 0.95, <po - 0, 
Oo = f , 6*0 = ^^0, 00 = - jc/o and ro = 
and ro = 20fi. In figure|5]we show some break- 



Fig. 3. The orbit with gravitomagnetic correction 
(red-hne) and without gravitomagnetic correction 
(black-Une). 

ing points of the orbital motion with gravito- 
magnetic corrections. 

3. Gravitational wave in the 
quadrupole approximation 

Now, considerin g th e orbital equations (see 
S. Capozziello et aT] (|2009|) ). we know 
tha direct signatures of gravitational radi- 
ation are its amplitude a nd its wave-form 
dC. W. Misner et al. I Il973h . In other words, 
the identification of a GW signal is strictly 
related to the accurate selection of the shape of 
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Fig. 4. The phiase portrait with (red-line) and with- 
out (blue-line) gravitomagnetic orbital correction. 
The example we are showing was obtained solving 
the system for the following parameters and initial 
conditions: yu ^ IMq, E = O.95,0o = 0, 9o = f , 
^0 = -^,00,00 = --^ro and i-Q = - j^ and ro = 20yu. 
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Fig. 5. We show the first four orbits in the phase 
plane: the red one is labelled I, the green is II , 
the black is III and the blue is IV. As it is possible 
to see, the orbits in the phase plane are not closed 
and they do not overlap the orbital closure point; 
we called these features 'breaking points Mn this dy- 
namical situation, a small perturbation can lead the 
system to a transition to chaos as prescribed by the 
Kolmogorov-Amold-Moser (KAM) theorem. 



wave-forms by interferometers or any possible 
detection tool. Such an achievement could give 
information on the nature of the GW source, 
on the propagating medium, and , in principle, 
on the gravitational theory producing such a 
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Fig. 6. In this figure it is shown the initial orbit with 
the initial(squares) and final(circles) points marked 
in black. 

radiation. It is well known that the amplitude 
of GWs can be evaluated by 

hj\t,R)^—^&' , (7) 

R being the distance between the source and 
the observer and {i,k} -1,2, where Qij is the 
quadrupole mass tensor 

a; = ^m„(3xi,4- V^). (8) 

a 

Here G is the Newton constant, r„ the mod- 
ulus of the vector radius of the a - th par- 
ticle and the sum running over all masses 
nia in the system. At this point we com- 
puted the amplitude components with grav- 
itoma g netic correc tions in geometrized units 
dS. Capozziello et al.i (l2009h ). We performed 
the numerical simulations in two cases: i) a NS 
of \AMq orbiting around a Super-MBH ( e.g. 
Sagittarius A* IO^Mq) ii) a BH of IOMq orbit- 
ing around a Super-MBH. We considered the 
reduced mass ji - J"'^^^ . Computations are 
performed with orbital radii measured in mass 
units. Initial distances are sampled to show or- 
bits from high eccentricity up to circularity 
(e = 1^"'"'^""" ). In Fig. |7]|9] we show respec- 
tively the field velocities of the orbits along the 
axes of maximum covariances, the total grav- 
itational emission waveform h and the gravi- 
tational waveform polarizations h+ and for 
a NS of \AMq. The waveform were com- 
puted for the Earth-distance from Sagittarius 
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Table 1. Data for GW for a NS of lAM^ orbiting around a Super-MBH 
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A (central Galactic Black Hole). We obtained 
the numerical examples solving the ODE sys- 
tem for the following parameters and initial 
conditions: fi 1.4Mo, E - O.95,0o = 0, 
00 = f,^o = 0, 00 = -j(/o and ro = 
and ro = (20 H- 1000)//. See also Fig. EH and 
Table 1. 
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Fig. 7. Plot of field velocities of the orbits. 



Finally we show in Fig. [TO] the plot of the 
estimated h GW-strain-amplitudes for the con- 
sidered binary sources at Galactic Center dis- 
tance. The blue-line is the foreseen LISA sensi- 
tivity (one year integration + white dwarf back- 
ground noise). The red diamonds ( 1 .4MO ) and 
the green circles ( 1OM0) are the h values for 
the systems we have studied. 




Fig. 8. Total gravitational emission waveform h for 
a neutron star of 1.4M„ 
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Fig. 9. The gravitational waveform polarizations 
h+ and for a neutron star (NS) of 1.4Mo 
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Fig. 10. Plot of estimated mean values of gravita- 
tional emission in terms of strain h for two binary 
sources from the galactic center with reduced mass 
ratio ji » IAMq (red diamonds) and fi « lOM© 
(green circles). The blue line is the foreseen LISA 
sensitivity curve. The waveforms were computed for 
the Earth-distance to Sagittarius A (central Galactic 
Black Hole). 
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4. Concluding Remarks 

The gravitomagnetic effect could give rise to 
interesting phenomena in tight binding systems 
such as binaries of evolved objects (neutron 
stars or black holes). The effects reveal partic- 
ularly interesting if v/c is in the range (10"' ^ 
10~^)c. They could be important for objects 
captured and falling toward extremely mas- 
sive black holes such as those at the Galactic 
Center. Gravitomagnetic orbital corrections, 
after long integration time, induce precession 
and nutation effects capable of affecting the 
stability basin of the orbits. The global struc- 
ture of such a basin is extremely sensitive to the 
initial radial velocities and angular velocities, 
the initial energy and masses which can de- 
termine possible transitions to chaotic behav- 
ior. In principle, GW emission could present 
signatures of gravitomagnetic corrections after 
suitable integration times in particular for the 
on going LISA space laser interferometric GW 
antenna . 
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